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Abstract 

Dynamical fluctuations in classical adiabatic processes are not considered by 
the conventional classical adiabatic theorem. In this work a general result 
is derived to describe the intrinsic dynamical fluctuations in classical adia- 
batic processes. Interesting implications of our general result are discussed 
via two subtopics, namely, an intriguing adiabatic geometric phase in a dy- 
namical model with an adiabatically moving fixed-point solution, and the 
possible "pollution" to Hannay's angle or to other adiabatic phase objects 
for adiabatic processes involving non-fixed-point solutions. 
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1. Introduction 

The adiabatic theorem is important in both classical and quantum me- 
chanics [H- It predicts a system's dynamical behavior subject to slowly 
varying system parameters. Although a general and mathematically rigor- 
ous proof of the adiabatic theorem is not obvious in both classical mechanics 
and quantum mechanics, the adiabatic theorem has been widely used. In- 
deed, it is always highly useful so long as there exist two drastically different 
time scales. The adiabatic theorem has also led to the discoveries of Berry 
phase [2] and the classical counterpart, i.e., Hannay's angle jsj. 

We focus on the classical adiabatic theorem (CAT), but as shown be- 
low, some of our results can be applied to quantum systems as well. Our 
interest here is not in a rigorous proof of the CAT, but in dynamical fluc- 
tuations around what is predicted by CAT. As discussed below, the possible 
consequences of the fluctuations neglected by the conventional CAT can be 
far reaching. The motivation of considering the fluctuations is based on a 
simple observation. That is, CAT, whose proof is based on an average over 
fast-varying variables, only reflects a mean dynamical behavior. As such 
fluctuations on top of a mean dynamical behavior should exist in classical 
adiabatic processes. Though fluctuations should be intuitively smaller in a 
slower adiabatic process, their effects are accumulated over a longer time 
scale and hence might not vanish even in the adiabatic limit. For instance, 
in a few early studies 0, [sl, 0|, including the study of "Hannay's angle of 
the world" [y, |7(], the actual total change in canonical variables may depend 
on the smoothness of the evolving adiabatic parameters. This abnormal 
behavior was shown to be connected with subtle fluctuations in the action 
variables from their average behavior predicted by CAT. Clearly then, a gen- 
eral description of the dynamical fluctuations in adiabatically evolving and 
classically integrable systems should be of importance. 

We shall present in this work a general result that describes the dynamical 
fluctuations inherent to classical adiabatic processes. Roughly speaking, it 
establishes an interesting connection between the actual rate of change of 
slowly varying system parameters and the actual classical orbits deformed 
from that predicted by CAT. To illustrate the usefulness of our general result, 
we design a simple dynamical model with an adiabatically moving fixed-point 
solution, from which an intriguing classical geometric phase can emerge. We 
then exploit our general result to discuss the "pollution" to Hannay's angle 
in classical adiabatic processes. A mean-field model that describes a two- 



2 



mode Bose- Einstein condensate (BEC) is also proposed to study fluctuation- 
induced "pollution" to adiabatic quantum evolution. 

To tackle with dynamical fluctuations, one may quickly think of an equa- 
tion describing the time dependence of the fluctuations around ideal adiabatic 
orbits. But this approach may not be fruitful because in principle, the time 
dependence of any canonical variables is already fully captured by classical 
canonical equations of motion. Instead, we are concerned with how fluctu- 
ations distort trajectories as compared with that predicted by CAT. In this 
sense, our approach is somewhat in a similar spirit as an early "multiple-time- 
scale-expansion" approach to corrections to classical adiabatic invariants in 
chaotic systems [8|. However, we focus on fluctuations associated with indi- 
vidual orbits in integrable systems, rather than fluctuations associated with 
an ensemble of chaotic trajectories in an energy shell. 

This paper is organized as follows. In Sec. II we derive a differential equa- 
tion describing the dynamical fluctuations in classical adiabatic processes. 
Some related details are also provided in Appendix. As an application, in 
Sec. Ill we study the case of an adiabatically moving flxed-point solution 
and show how an intriguing geometric angle may emerge in a simple toy 
model. Based on our general result. Sec. IV discusses why the "pollution" 
to Hannay's angle may exist and then proposes a physical system to study 
analogous fluctuation-induced pollution. We flnally give a brief summary in 
Sec. V. 

2. General Description of Dynamical Fluctuations in Classical Adi- 
abatic Processes 

Consider a classical integrable system with N degrees of freedom. Its 
Hamiltonian is given by H{p, q, R), where canonical variables p = {pi,P2, ■ ■ ■ , Pn) 
and q = {qi,q2, ■ ■ ■ , Qn) represent canonical momenta and coordinates, and 
R represents a collection of system parameters. Let F{I, q, R) be the gener- 
ating function that induces the R-dependent canonical transformation from 
(p, q) to the action-angle variables (1,0), where li = § Pidqi and = 
(6'i, 62, ... , On). 

To clearly present our derivation of a differential equation that describes 
dynamical fluctuations in classical adiabatic processes, this section is divided 
into four subsections representing the four steps in our derivation. First, after 
expressing classical equations of motion in the action-angle variables (I, 0), 
we deflne dynamical fluctuations on top of the idealized solution given by 
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CAT. Second, the time dependence of the canonical variables (p, q) is ex- 
pressed in terms of the dynamical fluctuations we define. Third, directly 
using the canonical equations of motion and the canonical transformation 
between the action-angle variables and the canonical variables, we reexpress 
the time dependence of the canonical variables in terms of the dynamical fluc- 
tuations as well as the action-angle variables along idealized classical orbits. 
Finally, by comparing results in the second and third steps a differential 
equation describing the dynamical fluctuations around idealized adiabatic 
orbits is obtained. 

2. 1 . Dynamical fluctuations 

In the (I, 0) representation an integrable Hamiltonian becomes R), 
which is independent of the angle variables 0. For time- varying R = R(t), 
the equations of motion for (I, 0) are given by |9|] 

dh _ _dW dR 

where a;i(I, R) = dl-i/dli is the angular frequency, and W is defined by 

W = VrF[I, q(I, 0, R), R] - p VRq(I, 0, R). (3) 

Note that Vr refers to the gradient in the parameter space under fixed (I, 0). 
If 



dR 



dt 



(4) 



is much smaller than c<;j|R|, one can take the average of Eq. ([T]) over the 
rapidly oscillating angle variables, yielding ^ ~ (W is a periodic function 
of 0). CAT hence identifies the action variables as adiabatic invariants, 
i.e., in adiabatic processes their values are fixed at I = (/i, /2, . . . , /at). For 
clarity, angle variables associated with this idealized solution are defined 
as = (6'i, 02, . . . ,0]\r). We also use p = p(T, 0, R) and q = q(I, 0, R) to 
describe the idealized solution in terms of the (old) set of canonical variables. 
With the action variables fixed at I, one may then solve Eq. ([2]) for a cyclic 
process from t = Otot = Tina straightforward manner. One may further 
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take the average of the ideahzed solution over all possible initial angle values 
to obtain Hannay's angle, which is the total mean angle change minus a 
dynamical angle. 

The above discussion does not represent a complete description of classical 
adiabatic processes. Clearly, Eq. tells us that ^ is not mathematically 
zero: it may possess fluctuations of the order 0(e) (i.e., to the first order of 
e). As such, in performing an averaging procedure as is done in CAT one 
neglects the dynamical correlation between and I. It is hence necessary 
to reconsider Eq. ([1]) in order to consider any possible real-orbit fluctuations 
on top of CAT. On a real orbit we assume we have /j = Jj + Sli, where we 
have used 6 to represent fluctuations from the behavior predicted by CAT. 
Equivalent to that, one can describe the same fluctuations from the idealized 
orbit in terms of 6qj and 6pj. 

There are now both idealized adiabatic orbits without considering fluctu- 
ations and true orbits with fluctuations: the geometry of an idealized orbit 
can be characterized by I = I and G [0, 2tt); and that of a true orbit with 
fluctuations is slightly deformed to 

I = I + 51, (5) 
= + 50, (6) 

where 61 and 50 are assumed to be at most of the order 0(e). By our 
definitions above, we have 

Here and in the following the summation convention for repeated indices 
is adopted. Equation ([7]) also defines two N x 2N matrices K and M, 
corresponding to the upper and lower halves of a Jacobi matrix. Note that 
throughout we use ^ to indicate ^ evaluated at x = x. 

° ax ox 

As will be seen below, it suffices to consider fluctuations of the first or- 
der of e because higher-order effects cannot be accumulated with time. We 
stress that the fluctuations are intrinsic: they are nonzero so long as e is 
not identically zero. In other words, fluctuations considered here exist in 
any classical adiabatic process and should not be thought of an effect arising 
from a too-large e. It should be also noted that in principle, all the dynamical 
information is contained in Eqs. ([T]) and ([2]). However, we are interested in 
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developing a framework to describe how fluctuations might behave along an 
idealized classical orbit. 



2.2. Canonical equations of motion in terms of fluctuations 

In terms of the fluctuations 6qj and 6pj, we next expand //(p,q, R) 
around H = if(p, q, R) to the order 0(e), yielding the following canoni- 
cal equations of motion for (q, p): 

dpi _ _dH _ d'^H _ 

dt dqi dqidpj ^ dqidqj ^ 

= ^.,(I,R^-^5,^-^^g,; 

dqi _ OH ^ ^'^^ Sp.^ 6q- 
dt dp I dpi&Pj ^ dpidqj ^ 

dq^ dm ^ dm ^ 

where we have used the following two canonical relations 

dij _ dpi _ 



dqi d9j' 
dij _ dqj 
dPi 86 ■ 



'0 



(9) 



Through Eq. ([8]) it is seen that the time dependence of the canonical variables 
(p, q) is connected to the dynamical fluctuations bq^ and to the first order 
of e. 



2.3. Time- dependence of canonical variables from action- angle variables 

The time evolution of the canonical variables (p, q) may be also directly 
obtained from the canonical transformation from the action-angle variables 
to (p, q) and from the equations of motion given by Eqs. ([T]) and ([2]). In 
particular, using 

dph _ dpi_dR dp^ 

dt ~ OR dt dij dt 39 j dt ^ ' 
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and the analogous expression for rewriting the derivatives in Eqs. (II]) and 
(|2]) at (I, 0) in terms of those at (I, 0), and neglecting all terms that are at 
least O(e^), one arrives at (see Appendix for details) 



dt 



dqi 
dt 



dp^dR _ dp^dW_ dR 
dR^ dlj dOj ' 'dt 
dW dR 



d6pi 
86, 



^ujj[l,R) 



09, 
dqi dR 
dR dt 

89 



= - + ujAVR) 



dlj dt 
dq^ dW 



dW 



dli 



09 
dR 



dR 

It 



35 qi 



39, ' 



XR) 



= - + ujAVR) 



3h 



Interestingly, due to the direct connection between (p, q) and (1,0), the 
full time dependence of (p, q) is connected with dynamical fluctuations in a 
highly nontrivial manner. In particular, the terms and in Eq. ffTTj) in- 

dicate that it is important to account for how dynamical fluctuations change 
with 9j. This is a crucial piece of information regarding the overall feature 
of the dynamical fluctuations. 

2.4- A differential equation describing dynamical fluctuations 

Both Eq dH]) and Eq. (fTTj) deal with the same time dependence of (p, q) 
and hence they should be consistent with each other. Comparing these two 
equations term by term, we arrive at the following equation, 

/ 5/i \ 

6h 



5p 
5q 



dSp 



^5^ 



(12) 



where 



(13) 



is a 2N X 2N matrix; 



is a 2N X 1 vector along each direction of R; and 




T-r _ ddj dh dOj dl2 /-, r\ 



is a 2N X matrix. Substituting Eq. ([7]) into Eq. (fT2|) . one finally obtains 
an equation for (5p and (5q only: 

V (nir - r) f ? Vs-^ = o. (16) 



For a given integrable Hamiltonian, except for those related to (5p and 5q 
and their derivatives, all the matrices contained in Eq. ( !T6|) are evaluated at 
an idealized orbit and hence can be explicitly obtained. 

Some remarks are in order. First, Eq. f lT6|) is not about evolving the fluc- 
tuations ((5q, (5p) at one moment to the next moment. Instead, it describes, 
when the system parameters reach the current configuration R with a small 
but nonzero rate the deviation of the overall shape of one true orbit from 
the idealized orbit without dynamical fluctuations, i.e., the overall deformed 
orbit in phase space. To our knowledge, this result is obtained for the first 
time here. This detailed description of the dynamical fluctuations can be very 
useful for both quantitative and qualitative considerations. The derivation 
here is somewhat lengthy because the physical meaning of (I, 0) in terms of 
(q, p) and hence the idealized orbit itself is changing as R varies. Second, 
consistent with our treatment to the first order of e, (5p, 5q) is seen to de- 
pend on If ^ were identically zero, then 5p = 5q = is one possible 
solution (If 5p 7^ and 5q 7^ is still the solution for ^ = 0, then this so- 
lution describes the relationship between two infinitely close orbits). Third, 
in the absence of the detailed information of ((5p, (5q) for at least one phase 
space location, Eq. ( !T6|) alone does not suffice to predict (5p, 5q) because of 
its differential form. As will be discussed later, this implies that in general, 
detailed information of the time-dependence of R, e.g., its smoothness, can 
be important for determining the dynamical fluctuations. Finally, because 
the linear Schrodinger equation and nonlinear Gross-Pitaeviskii (GP) equa- 
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tion have an exact canonical structure of Hamiltonian dynamics [10|, 111 
results here can be also relevant to quantum adiabatic processes. 

If we now consider the mean behavior of (5p, 5q) along an ideal orbit 
(denoted by (■)), then using the fact that (5p, 5q) are periodic functions of 
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0, we reduce Eq. (fTB]) to 

((n/.-r)(^,p)).(i:).^^o. (17, 

Because the matrices 11, K, T vary along the orbit, one may infer from 
Eq. ( !T7|1 the statistical correlations {{JIK — r)5p) and {(JIK — r)(5q), but 
the mean fluctuations (5p), (^q), or (51) remain unknown. 



3. Emergence of a geometric angle from an adiabatically moving 
fixed-point solution 

As a direct application of our central result in Eq. f lT6|) . here we focus on 
a rather simple case, where the solution to Hamilton's equation of motion is 
a fixed point in phase space if the system parameters are not changing. We 
denote the fixed-point solution as (p, q), which are of course functions of R. 
Consider now an adiabatic process in which R is changing slowly. Then the 
idealized orbit according to CAT is just one adiabatically moving fixed point. 
In addition, at this fixed point all functions of (p, q) are independent of © 
(otherwise they would be time-dependent), thus forcing their derivatives with 
respect to to vanish and making an averaging over [e.g., in Eq. (fTTll ] 
unnecessary. We therefore obtain 

n = 0, (18) 



(20) 



Using these results we have the following relation from Eq. (IT6 



\^) dt 



Note that the values of 9i at a fixed point can be taken as arbitrary. Hence 
the fluctuations obtained in Eq. (1201) do not have any interesting consequence 
for the evolution of 6',. Furthermore, since the K matrix vanishes at fixed 
points (where the action reaches its minimum), one would also arrive at 
51 = to the first order of e even though 5q 7^ and 5p 7^ 0. 

Consider then the coupling of this system with another degree of freedom, 
whose canonical coordinates are denoted by (J, 0). The total Hamiltonian is 
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assumed to be independent of 0, denoted iir*°*(p, q, J). Because J is a strict 
constant of motion and can be regarded as a fixed system parameter for the 
motion of (p,q), the expression for 6p and 5q in Eq. f l20|) still applies to 
fixed points in the phase space of (p, q). To seek how fluctuations predicted 
by Eq. fl20|) may affect the motion in 0, let us now examine the angular 
frequency associated with 0, i.e., 



wj(p,q, J) = 



OH 



tot 



dJ 



(21) 



Clearly, the fluctuations Sp and 6q will lead to 

c I dujj{p,q,J) dujj{p,q,J) 

(5wj(p,q, J) = — 6p + — Sq. 



dp 



dq 



(22) 



This fluctuation in uj{p, q, J) induces an correction to the evolution of 0. 
Using Eq. fl2I]]) . one flnds an explicit expression for this correction as follows. 



'duj 
dp 

dp ' 



5p + 



dujj 
dq 



5q 



dt 



9^ )r-i 



dp 
an 



■ dR. 



(23) 



As seen from Eq. ( 123|) . obtained above no longer depends on T (so it will 
not vanish even in the e — )■ or T — )■ +oo limit). Rather, it depends on the 
geometry in the parameter space only. (f)^°^^ is hence identifled as a geometric 
angle that arises from the fluctuations in a classical adiabatic process. This is 
particularly interesting because here 51 = 0, i.e., even when the fluctuations 
in the original action variables are vanishing, there can still be a physical 
effect on another degree of freedom due to the dynamical fluctuations. 

To illustrate the result in Eq. f l23|) we have designed a simple toy model 
with two degrees of freedom in total. Speciflcally, the total Hamiltonian is 
given by 



H"^\pi,qi;J) = aJ+- 



ZL_ T 

X2 ^ 



(24) 



with R = {X > 0,y > 0), being a cyclic angular coordinate that forms a 
canonical pair with J, and a being a free parameter. For the {pi,qi) degree 
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of freedom, this system has a R-dependent fixed point 



qi = V^Y; (25) 
p, = /JX, 

where J represents a conserved value of the variable J. 

To calculate the fluctuation-induced geometric angle seen in the evolution 
of 0, note first 

fifftot 2 2 



and _ _ 

r)/,i T 9t). fl/,1 T On- 

(27) 



duj 2pi duj 2q^ 



One may also easily obtain that the matrix F here is just a 2 x 2 matrix, i.e., 

i 2X2 - I Q2H I - I 4J Q ) ' l^^J 



and 



where X and F are unit vectors along the X and Y coordinates. Finally, sub- 
stituting these intermediate results into Eq. ([23]), one finds the fluctuation- 
induced geometric angle 



9pi dq^ ) ^2x2 ( ^ 1 ■ rfR 



= i y (FrfX - XdY) = - jj dS. (30) 

as=c 

As seen from the above result, here the geometric angle induced by the fluc- 
tuations in the first degree of freedom may be interpreted as the flux of an 
effective "magnetic charge" uniformly distributed on the (X, Y) plane. The 
emergence of such a new classical geometric angle from our simple calcula- 
tions is hence intriguing. It should be emphasized that in obtaining 0'^°'''' 
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in Eq. (150]) . we did not seek new action-angle variables (/i,6'i) and (/2,^2) 
such that becomes a function of Ji and I2 only. Indeed it can be highly 
complicated in general to find such a new representation due to the coupling 
between the two degrees of freedom. This indicates that (/)'^°'''' here has a 
different meaning than Hannay's angle, because it represents a geometrical 
correction to the evolution, not to the evolution of the yet-to-be-found new 
angle variables 9i or 62- 

We also note that our result here is consistent with one of the found 
terms in the previous study of the so-called "nonlinear Berry phase" based 



on GP equation [12[. In particular, the GP equation considered in Ref. [12 
can be mapped to that of a classical Hamiltonian with two degrees of free- 
dom, with the nonlinear eigenstates mapped to classical fixed points (see also 
Sec. IV-B). Adopting our perspective here, the geometric phase contributed 
by deviations from nonlinear eigenstates as analyzed in Ref. 12| may be un- 
derstood as a classical geometric angle due to intrinsic fluctuations in classical 
adiabatic processes. Indeed, we have checked that if we apply Eq. f l20|) to 



the model considered in Ref. 12 , then we can obtain a fluctuation-induced 



geometric phase term that is identical with a Berry-phase correction term 



discovered in Ref. [12|. Note however, the focus of our perspective is on a 
general description of the important dynamical fluctuations in a broad class 
of classical adiabatic processes. In our fully classical considerations here, a 
totally classical geometry angle is shown to arise in a second degree of free- 
dom that is coupled with the first degree of freedom (with one adiabatically 
moving fixed point solution); whereas in Ref. jl2|, the emphasis was placed 
on a quantum adiabatic evolution context and the main concern is with the 
sum of one familiar Berry phase and a fluctuation- induced geometric phase 
as a correction. 



4. Discussion 

4.1. Pollution to Hannay's angle 

As mentioned above, in some early studies about Hannay's angle in some 
Hamiltonian systems [3, 0, 0, 0|, it was numerically found that during an 
adiabatic process the total angle change minus the dynamical angle may 
not be Hannay's angle. This subtle behavior was connected with dynamical 
fluctuations in classical adiabatic processes. Here we exploit our general 
result of Eq. f lT^ to shed more light on possible pollution to Hannay's angle. 
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According to Eq. ([2]) and CAT, the total change in angle variables in a 
cyclic adiabatic process is given by 



On the right hand side of Eq. fl3Tl) . the first term is often called the dynam- 
ical angle, and the second term gives Hannay's angle (upon an average over 
initial angle variables). We have also used the notation 0'^^'^^^ to emphasize 
that it is for idealized cases without considering any dynamical fluctuations. 
Indeed, the angular frequency Ui in Eq. (13T1) is naively assumed to be the 
one determined by the idealized and constant action I. 

However, as suggested by Eq. ([2]), fluctuations in the action variables 51 

can then correct the angular frequency from a;i(I, R) to a;j(I, R) + ^'^gj^-' ■ 51. 
In terms of the canonical variables (p,q), fluctuations in p and q will lead 
to fluctuations in the angular frequency 

. .T 5W(I,R) e 5W(I,R) , 

5a;(I, R) = — ^ ■ 5p + — ^ ■ 5q. (32) 

For this reason, the dynamical angle obtained by a time-integral of the ide- 
alized frequency Ci;j(I, R), [see Eq. ( 13T1) ] should be re-examined with care. In 
terms of (5p and (5q, the real change in the angular variables should be given 
by 



^^real^y^ = A^^'^*^"'(T) + / 5uJi{l,K) dt 





Jo <9p 

+ r^.iq,, (33) 







Because 5p, (5q and hence 5uj are of the same order with e = \d'R,/dt\, just like 
the above fixed-point solution case, the term jj" 5uji dt may not be negligible 
as it accumulates the fluctuations 5uji{l, R) over an entire adiabatic process. 
So the term 5uji dt should not be neglected without a clear understanding 
of the dynamics. At this point it is also clearer why we only consider 5p 
and 5q to the first order of e: including higher-order terms are unnecessary 
because they will vanish in the e — )■ limit. 
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The correction term 6uJi dt can hence give the difference between two 
objects: the standard Hannay's angle, and a numerical calculation of a ge- 
ometric angle based on the expression of (A^^'^^^ - /q oJiil.Yi) dt). Unfor- 
tunately, unless for special fixed-point solution cases analyzed above, we in 
general cannot determine the fluctuations 5p and 5q from the differential 
equation in Eq. flTB]) . In particular, 5p and 5q can only be determined if 
we have information about them for at least one given (as the input). 
Therefore, without some detailed information of an adiabatic process, e.g., 
the detailed dependence of adiabatic parameter R(t) on time, information 
about 5p and 5q is not available in general. 

To see more clearly, let us discretize the adiabatic process by dividing one 
adiabatic process into many time intervals ti,t2,---, during each of which 
R = Rj, followed by a jump onto the next value Rj+i after the temporal in- 
terval tj (different time intervals and different choices for Rj define different 
adiabatic processes with different details). Note that even for a continu- 
ous adiabatic process, this discretized version is rather typical in numerical 
simulations (as the discretized time steps decrease, the simulated dynam- 
ics approaches a continuous process). Now for each point Rj, we may use 
Eq. (fT6|) to describe the dynamical fluctuations, but Eq. ( [T6|) is dependent on 
Rj. For a particular segment where R = Rj, the angle variable changes 
rapidly. Obviously, different timing for the next jump will result in different 
initial values of for next segment R = Rj+i, leading to another initial 
condition for the differential equation (ITB]) associated with R = Rj+i. This 
process then continues. According to Eq. fl32|) . 5uj and thus the correction 
term 6oOi dt will then depend on great details of a particular adiabatic 

process. It is for this reason that the correction term 6uji dt is identified 
as "pollution" to Hannay's angle, with the latter independent of how an adi- 
abatic process is implemented. Analysis here also makes it clearer that the 
fixed-point solution case in Sec. Ill is special because a definite prediction 
about fluctuations can be made therein. 

It is also worth noting that, according to Eq. fl32|) . the pollution vanishes 
if the angular frequency Ui does not depend on the action I. This is the case 
in a linear system such ClS db harmonic oscillator. 

4-2. "Pollution" to adiabatic phase evolution in a two-mode BEC model 

Finally, we propose to use a two-mode GP equation to study pollution 
to a geometric phase associated with quantum adiabatic cycles, thus making 
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a connection between our theoretical considerations here and a reachable 
experimental context. In particular, there are a number of possibilities to 
experimentally realize a two-mode BEC. For example, one may consider a 
BEC in a double- well potential, or a BEC in an optical lattice occupying two 



bands [13|. On the mean-field level, a two-mode BEC can be described by 



the following CP equation [h = 1) 



1 /7 + c(|6p-|ap) A \fa 

2 I A -7-c(|6|2- lap) J I b 



(34) 



where 7 denotes an energy bias between the two modes, |ap and |6p (with 
|ap + |6p = 1) represent occupation probabilities of the two modes, c gives the 
self-interaction strength, and A denotes the coupling between the two modes. 
We can consider, for example, the two parameters 7 and A to implement an 
adiabatic cyclic process. 

The dynamics described by the above CP equation can be translated into 
Hamiltonian dynamics. In particular, let p = (pa ~ 'Pbi Q = Ic^P? Q- = |a|e*''^'', 
b = |6|e*'^'', then apart from an overall phase parameter (pfy, Eq. flM|) leads to 



dp _ OH 
dt dq ' 

dq _ dH 

dt dp 



(35) 



where 



It is also straightforward to find that the evolution of (pb obeys 



H = Av/^(r^+|(2g-l)-|(2g-l)2. 



(36) 



where 



A = -^(2g-l)^ (37) 



It is seen that the evolution of the overall phase 0b is determined by, but 
will not have a back action on, the classical trajectories determined by H 
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in Eq. fl36p . In this sense, the 0^ parameter plays a similar role as the (j) 
parameter in Sec. III. 

It is now clear that our general result of dynamical fluctuations in classical 
adiabatic processes can be directly relevant to understanding the adiabatic 
evolution of a two-mode BEC system. If the adiabatic process starts from 
a stationary state of the GP equation, then the dynamics is just about an 
adiabatically evolving fixed-point solution of the Hamiltonian in Eq. fl5^ . 
As shown earlier (see also Ref. [HI), in this definite prediction can be 

made about how accumulation of dynamical fluctuations can eventually lead 
to a geometry-like correction to (pb- Consider now a superposition state of 
two stationary states of the above two-mode GP equation as the initial state 
of an adiabatic process. This case then corresponds to a classical adiabatic 
process with non- fixed-point solutions. As indicated by Eq. f p6|) . dynamical 
fluctuations can now affect the evolution of the adiabatically evolving phase 

in an unpredictable way if we do not know the details of the adiabatic pro- 
cess. Pollution to the quantum phase 0b hence emerges. Interestingly, in the 
same context, how (pb may develop an adiabatic geometric phase f or g eneral 



superpositions of stationary states was already considered in Ref. [IJ] with- 
out considering dynamical fluctuations. It is hence of interest to numerically 
or even experimentally examine the actual pollution due to the accumulation 
of dynamical fluctuation effects in such type of quantum adiabatic processes. 



5. Summary 

To summarize, we have obtained a general description of the intrinsic 
dynamical fluctuations in classical adiabatic processes associated with inte- 
grable systems. These fluctuations are typically neglected by the conven- 
tional classical adiabatic theorem. The dynamical fluctuations are described 
in this work in terms of deviations from idealized adiabatic trajectories. As 
an application, we have shown how a new kind of classical geometric phase 
may emerge using an explicit example with an adiabatically evolving fixed- 
point solution. We then discussed the origin of the pollution to Hannay's 
angle and proposed to use a two-mode BEC system to further study possi- 
ble fluctuation-induced pollution to one type of quantum adiabatic evolution 
described on a mean-field level. 
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Appendix A. On the derivation of Eq. 



Here we present some necessary details in deriving Eq. ( ITTj) . The sum- 
mation convention by using repeated indices is also adopted here. First, by 
definition we have Pi =Pi + Spi, and hence 



dpi 
OR 



dpi d6pi 
dR OR 



(A.l) 



As a second step, let us expand the expressions ^ and ^ around ^ and 



aii 



to the first order of 51 and 50, leading to 



dij 9/, fc^^. 



dijdik 



and 



dpi 



dPi . d'^Pi XT , J^^Pixa 
+ „^ oh + - - Odk- 



09, dO.dh 



To proceed further we shall use the obvious two relations 



dpi 
lit 

u;,(/,R) 



dpi dR dpi dIj dpi dOj 
dR~dt ^ WjHi ^ 

u,{lR) + ^5h. 
oh 



d9j dt ' 



(A.2) 

(A.3) 

(A.4) 
(A.5) 



Substituting Eqs. (1), (2), (Al), (A2), and (A3) into Eq. (A4), we find 

dpidR dpidW dR 



dpi 
dt 



dR dt 



86, 



dW dR 



dt 



+ 



= - + ujJl,R) + ^6I, 

_dlj dt ' ' dh 

54 + -|!^5^,)a;^.(I,R), 



(A.6) 
(A.7) 
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with all the terms of the second order of ^ or higher neglected. Note that 
consistent with our final result, we have assumed that {6p, Sq) or [61, 50) 
are of the first order of e = 

As the last step we use the relation 

Plugging this simple relation into Eq. (A6), we obtain the first equality in 
Eq. ( ITTj) . The second equality in Eq. (ITTl) can be obtained in the same 
manner. 
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